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This paper calculates the amount of baryon/radiation isocurvature fluctuation produced through 
the decay of a curvaton field, ft is shown in particular that if curvaton decay preserves baryon number 
and the curvaton dominates the energy density at the time of decay, the initial curvaton/radiation 
isocurvature mode is entirely transfered into a baryon/radiation isocurvature mode. This situation is 
opposite to that previously studied in three fluid models of curvaton decay; this difference is related 
to the conservation of the pre-existing baryon asymmetry and to the efficiency of the annihilation 
of all baryon/anti-baryon pairs produced in the decay. We study in detail the relevant cases in 
which the curvaton decay preserves or not baryon number and provide analytical and numerical 
calculations for each situation. 
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I. INTRODUCTION 

The curvaton scenario [l|, 0, H, 13, [H, B 0, B @ is a variant of the inflationary scenario in which the field driving 
the accelerated expansion (the inflaton field) is not necessarily that which produces all of the primordial fluctuations. 
Another field (the curvaton field), through its decay, can seed part of (or even all of) the cosmological perturbations. 
Therefore, in the most generic situation, these fluctuations originate from two different sources and the possibility of 
having isocurvature modes arises. 

The existence of isocurvature fluctuations would lead to distortions of the multipole moments of cosmic microwave 
background anisotropies, as compared to pure adiabatic modes. It is thus possible to constrain the fraction of 
isocurvature modes using high accuracy measurements, and present-day constraints show that the contribution of 
isocurvature modes is sub-dominant, at least if the isocurvature components are considered separately [lol. ITTI. H3, 
E3, EH El, E3, El, [3 US HH, H3|- Therefore, the study of the production of isocurvature fluctuations in the curvaton 
scenario, although interesting per se, can also help us to constrain the free parameters describing the model. 

The phenomenology of the curvaton scenario has been studied in the literature in a variety of cases, in particular 
in multi-fluid configurations (see for instance Refs. [!, H, 0, H, [H, H3, Hf| Ull)- The purpose of the present paper 
is to apply the formalism developed in Ref. [25| to the particular case of a net baryon/radiation isocurvature mode 
generated through curvaton decay. One peculiar feature that will emerge from the present study is the fact that 
the curvaton may induce a maximal isocurvature mode even if it dominates the energy density at the time of its 
decay, provided its decay preserves baryon number. This feature stands in sharp contrast with previous findings 
which showed that for curvaton decay into radiation and another fluid such as dark matter, the decay of a dominating 
curvaton would erase any pre-existing isocurvature mode. We study this case in detail and show that this particularity 
is related to the conservation of baryon number and to the efficient annihilation of all bb pairs produced in curvaton 
decay (throughout this paper, "b" stands for a generic baryon, and b for its antiparticle, not to be confused with 
bottom and anti-bottom quarks). 

This paper is organized as follows. In Sec. HH we describe the model and formulate the equations of motion at 
the background and perturbed levels. In Sec. IIII1 we numerically solve these equations in two cases, namely when 
the decay is symmetric in baryons and anti-baryons and when it is asymmetric meaning that the production of a 
net baryon number becomes possible. We show that these two cases correspond to very different phenomenologies. 
Finally, In Sec. IIV[ we discuss and compare our main results and present our conclusions. 
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II. DESCRIPTION OF THE MODEL 



We consider a model where four fluids are present: baryons (denoted in what follows with the subscript "b"), 
anti-baryons ("b"), radiation ("r") and the curvaton field ("cr"). At the fundamental level, the curvaton is a priori a 
massive scalar field but can effectively be treated as a pressureless fluid. One assumes that it can decay into radiation, 
baryons and anti-baryons. Each of these processes is controlled by a partial decay width denoted r CTr , and 
respectively. The curvaton decay occurs when the condition ~ H is met, where T a is the total decay width, namely 
T a = T al + r CT b + r,jb arLCl H the Hubble parameter. 

We do not discuss the phenomenology of curvaton to dark matter decay in the present paper. It is fair to assume 
that by curvaton decay, dark matter is effectively decoupled from radiation and baryons/anti-baryons. Even though 
the decay of curvaton may induce a dark matter - radiation isocurvature mode (see Ref. [25| for a detailed analysis), 
or a baryon - radiation isocurvature mode, both baryonic and dark matter sector will evolve independently. In this 
sense, the constraints obtained on dark matter or baryon isocurvature modes give complementary constraints on the 
physics of curvaton decay. 

The freeze-out of baryon/anti-baryon annihilations is controlled by the velocity averaged cross-section 

(o- b b«> - m~ 2 , (!) 

where m, = 135 MeV. This relation originates from the fact that, in the present context, the pion can be viewed as 
the gauge boson mediating the strong force. Freeze-out of bb annihilations occurs when r b |f = (cr^v) and/or 
Tb| f = rib|f ( (J bb v ) are °f the order of the expansion rate H which corresponds to a temperature ~ 20 MeV in the absence 
of curvaton decay (that is to say, assuming that radiation always dominates the energy content of the Universe). 

Big-Bang Nucleosynthesis (BBN ) pu ts rather extreme upper bounds on the amount of energy density injected at 
temperatures T < 1 MeV (see Ref. [27| for a recent compilation). For all practical purposes, it suffices to impose that 
Td > 1 MeV to satisfy these constraints. Furthermore, the late time decay of a scalar field at temperatures of order 
1 — 10 MeV is a fairly generic case in the framework of moduli cosmology. As is well known, such fields generically 
possess a very large energy density and a very small decay width T a ~ m^/m P1 , hence they decay after big-bang 
nucleosynthesis if their mass is of the order of the weak scale. Therefore, in order to reconcile the existence of such 
fields with the success of big-bang nucleosynthesis, one has two choices: either the energy density of these fields at 
the time of big-bang nucleosynthesis is very small or their mass is large, leading to early enough decay. The mass also 
cannot be arbitrarily large, otherwise one has to face a hierarchy problem, hence the generic decay temperature is 
1 — 10 MeV. Supersymmetric models with anomaly mediated supersymmetry breaking provide an explicit realization 
of particle physics model building in which the masses of moduli is of the order of m a ~ 10 — 100 TeV, which leads 
to decay temperatures of the moduli/curvaton Td ~ 1 — 10 MeV [28| . 

The above motivates the present study of the phenomenology of curvaton decay at temperatures of order 1 — 10 MeV. 
Out of simplicity, we keep this temperature fixed to a value T<j = 5.9 MeV in our numerical analysis, which corresponds 
to a total decay width T a — 1.6 x 10~ 20 MeV. We will argue that the results obtained remain unchanged if the decay 
temperature is higher, in particular if Td > 20 MeV. 

At the background level, following the approach of Ref. 25], the above situation can be modelled by the following 
set of equations 
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r CT b 




3 (a hb -v) 


TOp! H 


H 


8tt 


m h 






3 (a hi v) 


m\ x H 


H 


8tt 


m 5 



= n r n 6 + t^n a - " X " bb a " / '^ P1 — (n b n 5 - f^ q ) - (3) 



diV 

dfib 

dN 

H = (fir - 1) fir + ^ + 2 3Kb 8 ^ Wp1 ^ (fibfib " fib q fi C b q ) > ( 4 ) 

|k = n r o CT - ^ , (5) 

< W (6) 



dN 2 V 3 

Let us describe these equations in more detail. As usual, the parameters f2( Q ) are defined as the ratio of the energy 

density of the fluid a to the critical energy density, ^(a) = P(a)/Pcr- The time variable is the number of e-folds, 
N = In a, where a is the scale factor. The quantity f2 b q is defined by f2 b q = ra^n^ / p^, where n b q is the particle 
density at thermal equilibrium, expressed as: 

cq f m h T\ 3/2 ( TO b -/! b A ftT s 

n b = 9\ -7T- exp , (7) 



V 2?t J "V T 



3 



with a similar expression for rig . The quantity fib is the chemical potential of the baryons and one has /i b = — Mb- 
The temperature T can be expressed in terms of the variables of the previous system of equations as: 

^(^^y'V^. (8) 

\ 30 3mp[/ 

Note that the above description implicitly assumes that the curvaton decay products thermalize instantaneously. This 
assumption will be discussed at the end of Section UTTl One should already underline that the above ratios T^/Fa, 
r^b/TV and T^/T^ should be understood as characterizing the fraction of curvaton energy that eventually goes into 
thermalized "b", "b" and "r", rather than the branching ratios associated with curvaton decay channels. 

For the sake of simplicity, we ignore any temperature dependence of the function g* and we take — 10.75. If 
we compare with the equations of motion established in Ref. (25| in the case where the curvaton can decay into dark 
matter \ (rather than baryons and anti-baryons), the only difference is that terms like fl 2 , or ftT, cq are replaced 
by Ob^f, and f£ b q f2g q . Notice that, as a consequence, the evolution of the system does not depend on the chemical 
potential which cancels out, thanks to the fact that /ib = — Mb- Finally, there is a factor 2 in front of the last term in 
Eq. (U). This factor originates from the requirement that the total energy density be conserved. 

Let us also discuss how the initial conditions are chosen. Initially, we start with thermal equilibrium and some 
baryons/anti-baryons asymmetry. This implies that 



n b fi b = fi b eq) ft b eq) , n h -n h = s. (9) 



These two relations lead to 




Therefore, if the initial values of f2 b cq ^ and 5 are known, then one can deduce the initial values of fib and ilg. The 
quantities f2 b eq ' and S can be expressed as 

o(cg) _ 9 8tt mj 3/2 -x r _ 87rC(3)gmfc b 
b (2tt) 3 / 2 m 2 m|i ' 3H 2 m 2 pl n 2 x 3 ' { ' 

where x = m^/T and where the quantity £b is defined by 

e h ^^±. (12) 

n 1 

The present-day value of the baryon asymmetry is eb — 5.4 x 10~ 10 [29]. The initial value of t\, well before the 
freeze-out of bb annihilations must therefore be tuned in order to reproduce the final value after curvaton decay 
and entropy transfer from e + e~~ to the photons. Curvaton decay may dilute any pre-existing asymmetry through 
entropy production or even produce net baryon number if the curvaton decay process violates baryon number. In 
all our calculations presented further below, we have tuned this initial asymmetry in order to match the observed 
present-day value. 

In order to establish Eq. (TTJ]), we have used the fact that the number of photons is given by n 7 = ((3)gT 3 /n 2 . It is 
important to notice that the difference — is normalized with respect to the photon energy density (or number) 
and not to the total radiation energy density. In these formulas, x — x- m i ~ 10 (for instance) and ~ 0.9GeV are 
known (or chosen). Moreover, the Hubble parameter and f2 r are related through Eq. ©. Then, using the fact that 
the space- like sections are flat, i.e. f2 CT + Q r + ilb + f^b = L an d considering the (initial value of) fi^mi as a free 
quantity, one can derive the following expression 
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Therefore, for a given value of Slo-jni, i/mi can be computed and the other quantities f2 r ,ini, r2b,ini and f2 b ini , simply 
follow from the above equations. 

Let us now consider the perturbations. In order to establish the gauge-invariant equations of motion, we follow the 



method of Ref. 25j . It consists in formulating the equations in a covariant way in order to be able to perturb them 



consistently. One can write 

V /i T /i l/(ce ) = Q„( a ) + y„( a ) (14) 
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where = YT^ v u v is the curvaton decay term, u v being the four velocity of a fundamental observer, and the term 
is a phenomenological description of the interaction term. It reads 



'bb< 



u p . (15) 



ra b 

Of course, a rigorous treatment of the problem would rely on the full Boltzmann equation but this phenomenological 
description will be sufficient for our purpose. In particular, one can check that Eq. (|14j) exactly reproduces the 
background equations (JU)-©. Moreover, it is straightforward to perturb Eq. (fT4"|) . This leads to the following system 
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where A Q = Ap a /p a is the gauge-invariant density contrast for the fluid a. The quantity A oq is defined by the 
following expression 

A^iQ + ^A... (21) 

With this definition, it is easy to see that we deal with a "closed" system of equations since x must be viewed as a 
function of H and il r , see Eq. ©. Let us now turn to the discussion of the solutions of the two systems of equations 
presented in this section. 



III. RESULTS 



The main parameters that govern the cosmological consequences of curvaton decay into radiation and baryon/anti- 
baryons are: (i) the time of decay of the curvaton, which is encoded in the total decay width IV, (ii) the respective 
branching ratios T a b/^a and T ah -/T a \ (in) the magnitude of the curvaton energy density at the time of decay, i.e. f2^ d 
when H = T rT . The main parameters are therefore the respective branching ratios and Jl< d . Note that the branching 
ratios are constrained by the measured baryon asymmetry e b ~ 5.4 x 10~ 10 . In particular the baryon asymmetry 
e b measured immediately after curvaton decay should be equal to 1.5 x 10 -9 , in order to obtain the measured value 
after the reheating of the photon fluid by electron/positron annihilations. According to whether r CTb — r CTb or not, 
two possibilities may arise. In the case of symmetric decay, meaning r CTb = r CTb , the baryon asymmetry is generated 
by some unspecified mechanism acting at a higher energy scale; it is simply diluted during curvaton decay by the 
extra entropy brought by the curvaton. In the case of asymmetric decay, the curvaton contributes to the net baryon 
asymmetry. We note that direct baryogenesis at a low temperature T ~ 10 MeV is very contrived; we will nevertheless 
study this case for the sake of completeness and discuss the robustness of the results for higher decay temperatures. 
These two scenarios indeed exhibit different consequences, as discussed in turn in the following. 

We will use the standard definition of the curvature perturbation in fluid (a) [3(| HH, H3| : 

C W --*-^--* + ?T =^, (22) 

P(a) 3 [1 + W( a )J 
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FIG. 1: Evolution of the background and perturbed quantities in the case in which the curvaton decays symmetrically into 
baryons and anti-baryons: r CTb = IVg = 5x lCP 4 r CT . The total decay width is r CT = 1.6 x 10 -20 MeV, and f2^ = 0.9 at a 
temperature T = 94 MeV. The top panel shows the evolution of Q r (solid red line), Q a (dotted green line), Ob (upper dashed 
dark blue line) and (lower dashed light blue line) . The middle panel shows the evolution of £ r (solid red line) and £b (dashed 
blue line). The bottom panel shows the evolution of the isocurvature transfer coefficient Sbr/S^}. The transfer of isocurvature 
perturbation is maximal as the radiation but not the baryon fluid inherits the curvaton perturbation. 

where, in order to express pt a )i we have not considered the interaction term. The corresponding definitions for the 
isocurvature modes read 

Sbr = 3 (Cb - Cr) , ^br = 3 (Cb " Cr) , S aI = 3 (&r - Cr) • (23) 

In particular, we will be interested in the transfer of the initial curvaton/radiation isocurvature perturbation into the 

final baryon/radiation isocurvature mode, as expressed by the ratio S^/S^h The quantity indexed with (f) [resp. (i)] 
is evaluated well after the decay (resp. well before). In the following, we also express quantities evaluated immediately 
before (resp. after) decay with the superscript <d (resp. >d ). 

A. Symmetric decay 

In this sub-section, we explore the phenomenology of models in which the curvaton decays symmetrically into 
baryons and anti-baryons, i.e. T a b — T a ^. We find that two situations may arise, according to whether the curvaton 
dominates the energy density at its decay, i.e. il< d ~ 1, or not. 

Consider first the case in which the curvaton dominates the energy density at decay, f2^ d ~ 1. We find that the 
transfer coefficient of the isocurvature mode is maximal, as exemplified for instance in Fig. [TJ The top panel of this 
figure shows the evolution of the background energy density in radiation (solid red line), curvaton (dotted green line), 
in baryons (top dashed dark blue line) and in anti-baryons (bottom dashed light blue line). The middle panel shows 
the evolution of the individual Q Q ) quantities, and the bottom panel the transfer of the isocurvature fluctuation. This 
latter shows clearly the emergence of a net baryon-radiation isocurvature fluctuation after the decay of the curvaton. 

The middle panel also reveals that in this case, the radiation fluid inherits the curvaton fluctuations (since C-y/C<? ~ * 1 
at T <C Td), while the baryon fluid remains unaffected. 

This result is quite different from a "standard" scenario of curvaton decay into radiation and dark matter, in which 
the domination of the curvaton at the time of decay ensures that only adiabatic modes subsist, as all fluids have 
inherited the same curvaton perturbations. This difference can be related to the annihilation of all bb pairs produced 
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by the curvaton, which effectively reduces to zero the net energy transfer of the curvaton to the baryon fluid. In order 
to put this statement on quantitative footing, it is useful to evaluate the ratio of the annihilation rates of baryons and 
anti-baryons to the expansion rate immediately after curvaton decay (H = ]7 CT ): 

= b \ bb / 16 d_ n > d 

T a \10MeVj b ' 

T _ = "b >d Kb^> „ 46 x 1Q i6 ( _g_y n>* . (24) 



WbbJ7 ^ 4 g x 1Q 16 ^ n>d 

r CT ~ ' V 10MeV , 

The first equation gives the ratio Tb of the annihilation rate of baryons to the expansion rate, while the second gives 
the corresponding ratio Tg of the annihilation rate of anti-baryons to the expansion rate. Considering Tb, the above 
formula shows that if fig exceeds ~ 10~ 16 , annihilations are effective. In the absence of curvaton, the freeze-out of 
bb annihilations occurs as the abundance of anti-baryons is reduced to below this threshold. In the presence of a 
curvaton however, the decay of this field will regenerate the annihilations provided the amount of curvaton produced 
anti-baryons is sufficient, i.e. 0> d - VL< d T ai /T a ~ T^/T a > 1CT 16 . Then all pairs of baryon/anti-baryon produced 
by curvaton decay will annihilate. Of course, if the branching ratio V a ^/T a < 1CP 16 , the regeneration of annihilations 
will not take place, but the curvaton will not exert any influence on the pre-existing baryon fluid either. 

In fact, the behaviors of the different quantities plotted in Fig. Q] can be understood in more detail along the 
following lines. Consider the variables associated to net baryon number, in particular Ob — Og. Its equation of motion 
reduces to: 

A (Ob - O b ) = O r (O b - O b ) . (25) 

This composite fluid is isolated, as neither annihilation nor curvaton decay violates baryon number. Therefore, the 
curvature perturbation associated to this composite fluid is conserved, as predicted in Ref. Q. Furthermore, one 
has O b <C Ob before curvaton decay and after annihilation freeze-out, so that this fluid of "net baryon number" 
approximately corresponds to the baryon fluid. As annihilations of bb pairs produced by curvaton decay is efficient, 
the above inequality remains valid after curvaton decay, hence "net baryon number" remains a good approximation 
for the baryon fluid. All in all, the above indicates that the curvature perturbation of the baryon fluid should remain 
conserved if curvaton decay preserves baryon number and if annihilations of bb pairs are efficient. 

This implies that the above results remain unchanged if the decay temperature Ta > 20 MeV. The theorem of 
Ref. stipulates that the isocurvature mode between two fluids sharing thermal equilibrium are erased on a small 
timcscale, unless there exists a conserved charge. In the present case, baryon number is conserved, or more precisely 
net baryon number does not couple to radiation, hence the above theorem does not apply. Consequently, once the 
isocurvature mode is produced, it remains conserved unless baryon number violating processes take place. In other 
words, one can extrapolate the above results to temperatures at least as high as the electroweak scale. Let us also 
remark that, if baryons are relativistic (at temperatures above the QCD scale), the above equations are slightly 
modified, but the above results remain unmodified. 

One a more formal level, one can follow the evolution of the different variables as follows. Neglecting O b in front 
of Ob in Eq. (|25[) above indicates that O b scales as a when O r ~ 1 (i.e. after curvaton decay), while Ob remains 
approximately constant when O r 1. These trends are observed in Fig. |TJ 

The behavior of O b is less trivial to obtain (but its cosmological relevance is also much less) . One can approximate 
Eq. ([31) with the following, after curvaton decay: 

dOg _3ml 1 (a hh v)H 
dN 8tt TO b 

where the term O r Ob has been neglected as the annihilations are dominant. Using the fact that H cx a~ 2 and Ob oc a 
after curvaton decay, one derives the following late time value of O b : 

Of * 0>" e"^ d . (27) 

Since Tg d takes enormous values of order 10 9 , the annihilations regenerated by curvaton decay essentially erase all 
trace of anti-baryons and the corresponding plateau cannot be observed in Fig. [1] because it is too small. 

Let us now turn to the perturbations and assume that curvaton decay is instantaneous. If the curvaton dominates 
the energy density before decay, and transfers its energy to radiation, then: 

(P - CP ■ (28) 
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FIG. 2: Same as figure Q] except that fi^ = 0.01, which corresponds to O. 



><d 



0.09. Essentially no isocurvature fluctuation is 



produced as neither the baryon nor the radiation fluctuations have been affected by curvaton decay. 



This relation can be obtained through standard methods and corresponds to the conservation of the total curvature 
perturbation throughout curvaton decay. Similarly, one can build the variable associated to the perturbation of net 
baryon number, fibAb — fibA^,, which for all practical purposes, can be approximated by fibAb. The equation of 
motion for this quantity reads: 



A (f] bAb - fi b A b ) = 3|| - 



fi r (fi b A b - fi b A b ) 



(29) 



Since <f> is conserved both before and after curvaton decay, the first term on the r.h.s. can be neglected, and fi b Ab — 
fi b A b is approximately conserved when fi r <C 1. Approximating fi b A b — fi b Ab with fi b A b , this implies that A b 
is approximately conserved, since fi b is constant in this case (see before) and, hence, that £b is also conserved. At 
late times, after curvaton decay, fi r ~ 1 implies that fibAb scales as a, hence that A b (and therefore £b) is again 
approximately constant because fi b oc a. One thus finds that: 



At) 

^b 



Sb 



(30) 



As mentioned above, this property can be traced back to the fact that net baryon number behaves in the present case 
as an isolated fluid, hence its curvature perturbation is a conserved quantity. Finally, one derives from Eq. (|28[l and 
(|3"0"|) above the transfer of isocurvature perturbation: 



(f) 



-5« (fi< d ^ 1) 



(31) 



These results match the numerical evolution observed in Fig. [T] 

Obviously, the above discussion suggests that S^} — > as fi< d — > since the net baryon number must remain 
unaffected, while a decreasing curvaton energy density at the time of decay implies that a lesser amount of radiation 
is produced during the decay. In more detail, one should obtain 



S, 



(f) 

br 



"fi< d ^ 



since 



CP * (i-fi< d )C (i) +fi< d C (i) 



This trend is confirmed in Fig. [5] which provides an example with fi^ = 0.01, corresponding to fi< d 
(Td — 5.9 MeV). The final transfer coefficient is of order — fi< d as expected. 



(32) 
(33) 

0.09 at decay 



8 



B. Asymmetric decay 



If the inflaton can decay asymmetrically, T a b ^ T a ^, the phenomenology is different, as all "b" and "b" produced by 
the curvaton will not be able to annihilate with each other. In particular, the production of net baryon number during 
curvaton decay comes with the transfer of the curvaton perturbations to the baryon fluid. As already mentioned, 
known models of baryogenesis produce baryon number at a much higher scale than 1 — 10 MeV. We nevertheless discuss 
this asymmetric case for the sake of completeness and because it provides useful insights into curvaton cosmology. 
Moreover, as we have argued in the previous section, the present results can be extrapolated to a higher decay 
temperature, possibly as high as the electroweak scale. 

In the present case, one may expect cosmological consequences opposite to those found in the case of symmetric 
decay: if the curvaton dominates the energy density of the Universe shortly before decaying, and produces during its 
decay most of the baryon number, both baryon and radiation fluid will inherit its perturbations, hence there should 
be no final baryon/radiation isocurvature mode. On the contrary, if the curvaton energy density is small compared 
to the radiation energy density shortly before decay, but the curvaton still produces most of the baryon number, a 
maximal isocurvature mode between baryon and radiation should be produced. 

These trends are confirmed by the numerical computations, as shown in Figs. [3] and 2J The first figure, Fig. [31 
corresponds to the same value of 0,$ as in Fig. [TJ but with an asymmetric decay width (T^ — T^/T^ = 2 x 10~ 8 . 
In what follows, we will use the short-hand notation: 

AB bb = T ^- T ^ . (34) 

Assuming that the initial baryon asymmetry vanishes and that curvaton decay is instantaneous, one can obtain an 
order of magnitude of the decay asymmetry needed to reach the observed value of 6b as follows: 

n>* - n>« ~ AB hh n<* , (35) 

which implies: 

eb - 7.3 x 10- ( 1Q _ 2 ] o CT MeV ) V2 AB h , fi<* . (36) 

Numerical calculations differ from this simple estimate by a factor of order unity. 

In order to understand these results, it is instructive to express the time evolution of the baryon asymmetry using 
the system of Eqs. ((2]) , (J3j) , (J4j) and (J6j) . The baryon asymmetry can indeed be written as: 

45 \ ;/t ( mp.y'7 fl-V " SV^s „ (H_ V" SV^s 



Hence the time evolution of the baryon asymmetry is governed by the following equation: 

i de b ^ &<j r gb - r CTb _ 3 fig- r gr 

e b diV ~ n h - n s H 4 ft r H 



(38) 



In order to obtain the above equation, we have neglected the baryon/anti-baryon annihilation term in the equation 
for f2 r [Eq. Q], which is justified insofar as the amount of radiation produced in baryon/anti-baryon annihilations is 
negligible at or after freeze-out. 

The above equation is interesting because it shows how the baryon number can be modified: either through baryon 
number violating curvaton decay (first term on the r.h.s), or through dilution due to entropy production (second term 
on the r.h.s). It also provides an estimate of the conditions under which the initial curvaton/radiation isocurvature 
mode is efficiently transfered to the baryon/radiation mixture. Such an efficient transfer can indeed be achieved if 
|Aeb/eb| ~ 1 at curvaton decay, without significant production of radiation by the curvaton. The latter condition 
amounts to negligible entropy production, or what is equivalent, to assuming that the second term on the r.h.s of 
Eq. ([55)) is negligible compared to unity. The former condition then implies that the first term on the r.h.s of Eq. ([3"5[) 
is larger than unity. All in all, efficient transfer of the isocurvature mode occurs if: 

AB bS n<<> > n<« - n<« , < n<«« . (39) 

It is interesting to remark that this situation is very similar to that encountered for curvaton decay in a three-fluid 
model incorporating radiation and dark matter. Borrowing from the method of Refs. [23l.[25|. it is possible to express 
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FIG. 3: Same as figure [T] except that curvaton decay now violates baryon number, with r CT b — = 1.7 x 10 8 r CT . The 
final baryon number matches the observed value, for an initial asymmetry tb = 10~ 13 . Other quantities remain unchanged, 
in particular f2^ = 0.9 and r CT = 1.6 x 10~ 20 MeV. Essentially no baryon/radiation isocurvature fluctuation results, since the 
baryon and the radiation fluctuations have been similarly affected by curvaton decay. 
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FIG. 4: Same as figure[3]for baryon violating curvaton decay, except that r^ 1 ' = 0.01, which corresponds to fla d = 0.09. The 
baryon violating decay width is such that the final baryon number produced matches the observed value; this corresponds to 
A£> b g = 1.3 x 10 -7 for an initial asymmetry eb = 10 -13 . A large isocurvature fluctuation is produced as the baryon (but not 
the radiation) fluctuations have been affected by curvaton decay. 
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the final baryon/radiation isocurvature fluctuation in terms of the initial curvaton and radiation curvature modes, as 
follows. One first constructs a composite fluid that has the property of being isolated, with energy density: 



-'comp 



= Ph 



Ph 



AB bbP<? 



(40) 



Notice that this construction is possible because each component of the composite fluid is pressureless. Its curvature 
perturbation, which is conserved by construction, is: 



Co 



-Cb 



-Cb 



"Cor 



(41) 



Then, assuming that curvaton decay is instantaneous, one can match the value of Ccomp after decay to that before 
decay, which gives: 



sb 



^ comp 



^comp 



of bb pairs after curvaton decay. Although the quantity Ccomp * s evaluated here immediately before decay, it can be 



In order to obtain the first equality, we have used the fact that Q 
of bb pairs after curvaton decay. Although the qi 
evaluated at any initial time, since it is conserved 

The radiation perturbation is given by Eq. 
written as: 

c(f) 
^br 



b 



"b 



(42) 



as a result of the efficient annihilation 



33|) . hence the final baryon/radiation isocurvature perturbation can be 



AB hi a 



(i) 







(i) 



n 



(i) 



C(i) 



(43) 



where we used the fact that 5, 



00 



j(i) 



br 



br 



0. As expected, the isocurvature transfer vanishes as il< d — > (since this 

also implies — ► 0). When fl< d — > 1, one can see that the first term in the bracket on the r.h.s. of the above 
equation also tends to one, and therefore the transfer coefficient of the isocurvature mode also vanishes. The initial 
isocurvature fraction is transfered efficiently only if the conditions expressed in Eq. (|39[) are fulfilled. Note also that 
in the limit A_B bb — > 0, one recovers the result of Section lill Al presented in Eq. (|32l) . 

Finally, a last point is to be made concerning the assumption of instantaneous thermalization of the curvaton decay 
products. If the center of mass energy y/s ~ (EEth) 1 ^ 2 for an interaction between a high energy particle of energy 
E and a thermalized particle of energy E t h is well above the QCD scale, then the ratio of rates of thermalization 
processes to bb producing ones is of the order of (a Gm /a s ) 2 -C 1. It is even less if s/s is smaller than the QCD 
scale. Therefore the above approximation is not strictly speaking justified. However the neglect of these additional 
interactions would not modify our conclusions, for the following reason. 

The only effect that could modify our conclusions is if one fluid (either radiation or baryon) were "contaminated" 
by the other fluid (respectively baryon or radiation) through the interaction of high energy particles produced through 
curvaton decay with thermalized particles. One typical example is given by the transfer of energy from the photon 
to the baryon fluid through 7 + 7th — * b + b, where 7 stands for a high energy photon. However, net baryon number 
does not couple to radiation, hence transfers of energy between these two fluids cannot take place after curvaton decay 
(provided this latter occurs after any baryogenesis event). 

Hence all the conclusions remain unaffected by these processes that occur between curvaton decay and thermaliza- 
tion. It is important to stress, however, that T aT /Y a , T^b/Tu and Y a ^/Y a should not be interpreted strictly speaking 
as the branching ratios of curvaton decay into radiation, baryons or anti-baryons, but rather as the fraction of curvaton 
energy eventually transfered into these fluids after all thermalization processes have occured. 



IV. CONCLUSIONS 



In this section, we recap our main results. We have studied the production of isocurvature perturbations in 
the curvaton scenario where the curvaton field can decay into radiation, baryons and anti-baryons. Two different 
cases have been considered. The first one is the symmetric case in which the curvaton/baryon decay width equals the 
curvaton/anti-baryon one, i.e. curvaton decay preserves baryon number. We have found that if the curvaton dominates 
the energy density before decay, then a baryon/radiation isocurvature mode can be produced. In the opposite situation 
in which the curvaton contributes negligibly to the total energy density immediately before decaying, the isocurvature 
mode vanishes. This result is opposite to the standard prediction of the simplest curvaton scenario in which any 
pre-existing isocurvature mode is erased by curvaton decay if this latter dominates the energy density at the time of 
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decay. This difference can be traced back to the conservation of baryon number and to the annihilation of all bb pairs 
produced during curvaton decay. 

One noteworthy consequence of the above is to forbid the liberation of a significant amount of entropy by a late 
decaying scalar field at temperatures below any baryon violating processes, such as is often invoked for the dilution 
of unwanted relics. 

Another consequence of the above is that a baryon-radiation isocurvature mode Sbi cannot co-exist with a (WIMP) 
dark matter - radiation isocurvature mode S xr , since the conditions to produce these modes are opposite to one 
another. Since S x h = S XI — Sbr, the existence of a baryon-dark matter isocurvature mode appears generic in this 
case (unless fl< d is so small at the time of decay that the curvaton exerts essentially no influence on dark matter and 
baryon perturbations). 

The asymmetric decay presents a different phenomenology. Since the curvaton decay does not produce the same 
number of baryons and anti-baryons, the annihilations cannot suppress all the baryonic decay product and, as a 
consequence, when the curvaton dominates at decay, the isocurvature perturbations are erased. In this case, most 
or all of the baryon and radiation fluctuations indeed originate from the curvaton. If the curvaton contribution to 
the energy density is smaller than unity at the time of decay, then radiation cannot be affected substantially, while 
the baryon fluid may be strongly affected; this situation results in a large baryon/radiation isocurvature fluctuation. 
In some sense, this case appears similar to the case of curvaton to dark matter decay studied in Ref. [25| |. Contrary 
to the previous symmetric case, non vanishing S*b r and S X1 - can co-exist. We note however, that baryogenesis at low 
scales (below the electroweak phase transition) is rather contrived. 

On more general grounds, the study presented in this article exemplifies how scenarios where scalar fields can decay 
at late times can be constrained not only at the background level, as it is usually done, but also by investigating 
the consequences at the perturbed level. It is clear that, if this type of information is taken into account, one can 
hope to improve our understanding of the feasibility of such theories. We hope to return to this question in future 
publications. 
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